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ABSTRACT 



Temporal and spectral characteristics of prompt emission of gamma-ray burst 
(GRB) pulses are the primary observations for constraining the energizing and 
emission mechanisms. In spite of very complex temporal behavior of the GRBs, 
several patterns have been discovered in how some spectral characteristics change 
during the decaying phase of individual, well defined long (> few seconds) pulses. 
In this paper we compare these observed signatures with those expected from a 
relativistically expanding, shock heated, and radiation emitting plasma shell. 
Within the internal shock model and assuming a short cooling time, we show 
that the angular dependence in arrival time from a spherical expanding shell can 
explain the general characteristics of some well defined long GRB pulses. This 
includes the pulse shape, with a fast rise and a slower decay, oc (1 +t/r) 2 , where 
t is a time constant, and the spectral evolution, which can be described by the 
hardness-intensity correlation (HIC), with the intensity being proportional to the 
square of the hardness measured by the value of the peak, e.g. E p of the vF v 
spectrum. A variation of the relevant time scales involved (the angular spreading 
and the dynamic) can explain the broad, observed dispersion of the HIC index. 
Reasonable estimates of physical parameters lead to situations where the HIC 
relation deviates from a pure power law; features that are indeed present in the 
observations. Depending on the relative values of the rise and decay times of 
the intrinsic light curve, the spectral/temporal behavior, as seen by an observer, 
will produce the hard-to-soft evolution and the so called tracking pulses. In our 
model the observed spectrum is a superposition of many intrinsic spectra arriving 
from different parts of the fireball shell with varying spectral shifts. Therefore, 
it will be broader than the emitted spectrum and its spectral parameters could 
have complex relations with the intrinsic ones. Furthermore, we show that the 
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softening of the low-energy power-law index, that has been observed in some 
pulses, can be explained by geometric effects and does not need to be an intrinsic 
behavior. 

Subject headings: gamma rays: bursts — gamma rays: theory — relativity 

1. INTRODUCTION 

The mechanism underlying the prompt 7-radiation in gamma-ray bursts (GRBs) is still 
an unsolved puzzle. There is, however, a growing consensus about some aspects of it. The 
large energies and the short time scales involved require the 7-rays to be produced in a 
highly relativistic outflow, an expanding fireball. In the standard fireball model 7-rays arise 
from shocks internal to the outflow at a distance of R ~ 10 13 — I0 17 cm from the initial 
source. The episodic nature of the outflow causes inhomogeneities in the wind (or shells) 
to collide and thus creating the shocks. These tap the bulk kinetic energy and transform it 
into random energy of leptons which radiate. The dominant emission mechanisms are most 
probably non-thermal synchrotron (Tavani 1996; Lloyd & Petrosian 2000) and/or inverse 
Compton emission (Panaitescu & Meszaros 2000), but there have been other suggestions, 
for instance, thermal, saturated Comptonization (Liang 1997). 

The fundamental process of a burst is thus an individual shock episode which gives rise 
to a pulse in the 7-ray light curve. Superposition of many such pulses create the observed 
diversity and complexity of light curves (Fishman et al. 1994). The spectral and temporal 
characteristics of these pulses hold the key to the understanding of the prompt radiation of 
GRBs. However, there is no consensus on what effects lie behind the observed pulse shapes 
and their temporal and spectral evolution. 

The overall spectra of most GRBs can be described by a simple broken power law with 
a low and a high energy index, say a and /3, and a break energy E^. Often a > —2 and 
(3 < —2, so that the vF v or EFe spectrum peaks at a photon energy E p ~ E^. The total 
light curves of GRBs, on the other hand, are very diverse and not readily describable by a 
simple formula. Nevertheless, many attempts have been made to decompose the complex 
light curves into pulses and analyze their characteristics (Norris et al. 1996; Lee et al. 
2000a,b). No simple patterns have emerged from these studies of the population as a whole. 
However, some relations have emerged from investigations of GRBs with simple light curves; 
those described by a single pulse or a few, well separated pulses (Kargatis et al. 1995; Ryde 
& Svensson 2000, 2002; Borgonovo & Ryde 2001) (hereafter BR01). The pulse shapes and 
evolution of spectra seem to obey some simple relations. Motivated by these results, in this 
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paper we explore possible explanations for these behaviors. 

Several different possibilities exist. The simplest scenario is to assume an impulsive 
heating of the leptons and a subsequent cooling and emission. The rise phase of the pulse 
is attributed to the energizing of the shell which we will refer to as the dynamic time and 
the decay phase reflects the cooling and its time scale. The instantaneous spectrum reflects 
the cooling of the lepton distribution. The primary problem with this interpretation is 
that, in general, the cooling time for the relevant parameters is too short to explain the 
pulse durations and the resulting cooling spectra are in drastic disagreement with the above 
observed form (Gissellini et al. 2000). A more plausible model is one where the pulse 
duration is set by the dynamic time of say the shell crossing, which could be much larger 
than the microscopic acceleration and/or emission-cooling times. In this case there is a 
continuous acceleration of particles during shell crossings; the acceleration and the cooling 
occur in situ and simultaneously and give rise to the observed behavior. The pulse shape 
then is a reflection of the energizing mechanism of the electrons. A third possibility is that 
the above picture operates only during the rise phase of the pulse and that the decay shape 
is due to geometric and relativistic effects in an outflow with a Lorentz factor of T ^ 100. 
The curvature of the fireball shell will make radiation, emitted off the line of sight (LOS, for 
short) delayed and affected by a varying relativistic Doppler boost, due to the different light 
paths the photons have to travel. 

The aim of this paper is to investigate to what extent, and how, the last model affects 
the observed light curve and spectral evolution during the individual pulses; in particular 
to determine whether the resultant behavior can explain the observed relations found for 
simple pulses mentioned above. We want to point out that the discussion of the observable 
signatures due to the curvature effect is independent of the process underlying the intrinsic 
pulses of radiation. Colliding shells and internal shocks are one example. However, other 
possibilities exist, for instance, as Lyutikov & Blackman (2001) pointed out, if the outflow 
is Poynting-flux dominated, the intrinsic radiation could be caused by a non-linear break- 
down of large-amplitude electromagnetic waves at a distance of approximately 10 14 cm from 
the progenitor. Furthermore, we emphasize that the description is for individual emission 
episodes, i.e. single pulses and one must bear in mind the possibility that these actually 
could consist of several heavily overlapping pulses; see further discussion in BR01 and Norris 
et al. (1996). The observations relevant to our discussion will be described in §2 and the 
appropriate time scales in §3. In §4 we derive the spectral and temporal structure expected 
in a simplified version of the proposed model. A more realistic model including both the 
dynamic and curvature effects is discussed in §5. Some other complications and caveats are 
discussed in §6 and a brief summary and discussion of the conclusions are given in §7. 
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In the following, primed quantities are evaluated in the comoving frame at rest with 
the outflowing material in the shock front. The rest frame will denote the inertial frame at 
rest with the progenitor. The cosmo logical time dilation and spectral redshift, which are 
constant factors for individual pulses and bursts, will be ignored. 

2. OBSERVATIONAL SIGNATURES 

In this section we describe the observational signatures of GRBs with simple light curves 
that we wish to explain. Most GRBs exhibit complex light curves and only a small frac- 
tion have sufficiently long (say longer than few seconds) and smooth pulses to allow for 
detailed temporal and spectral investigations. Many studies have analyzed such small sam- 
ples of bursts and pulses and drawn some important conclusions, which give us clues to the 
underlying physical processes for the creation of these pulses. 

2.1. Hardness Intensity Correlation (HIC) 

One of the most important observational signatures is the correlation between the bolo- 
metric energy flux, F^ \(t) = J F(E, t) dE and the hardness of the spectrum during the time 
evolution of individual pulses, where, F(E,t) is the flux of the energy (not photon) spec- 
trum. This relation, referred to as the hardness- intensity correlation (or HIC for short), was 
first discovered by Golenetskii et al. (1983), using the temperature from thermal spectral 
fits as the measure of the hardness. Subsequent observations have shown that the thermal 
spectra do not provide a good fit to a majority of GRBs. More recently BR01, representing 
the 'hardness' of the spectrum by the peak photon energy E p of the vF v = EF E spectrum, 
find the simple power-law relation 

Fboi = F ho \fl(Ep/E P fl) v , (1) 

where E Pt0 and Fb i,o are some fudicial values of the peak photon energy and the bolometric 
energy flux, usually taken at the beginning of the decay phase. The time evolution of the 
spectrum is mostly from hard and intense to soft and dim. Figures la and 2 show examples of 
the HIC for pulse decays observed by the Burst and Transient Source Experiment (BATSE) 
on the Compton Gamma-ray Observatory (CGRO). More examples can be found by BR01. 
The distribution of the power-law index rj of long GRB pulses is somewhat broad. The 
original study by Golenetskii et al. (1983), using a thermal fit and kT for E p , found the 
power law index to vary between 1.5 — 1.7. Kargatis et al. (1995) studied 26 GRBs with 
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prominent pulses and found a power-law HIC behavior for 28 pulse decays in 15 of these 
bursts. Their distribution was centered on rj = 1.7. BR01 studied a sample of 82 GRB pulse 
decays and found them to be consistent with a power-law HIC in, at least, 57% of the cases 
and for these found rj = 2.0 ± 0.7 (see their Figure 3, depicting the sample of 47 pulses with 
good power-law HICs.) An example of a pulse in this sample with rj < 2 is shown in Figure 
2d. Other such examples can be found in Figures 4c, 6c, f in BR01. In the BROl-sample, 11 
cases have rj > 2.0 by two standard deviations or more. All of these differ somewhat from a 
perfect power-law with most having a concave shape. Three examples of these are given in 
Figure 2a, b, and c. Here, whenever possible, the HIC relation for the rise phase is included 
(see Fig. 3 for the time interval used in the different cases). In Figures 2b and 2d we observe 
a monotonic evolution of E p for both the rise and the decay phase, which is characteristic 
for the so-called hard-to-soft pulses, while Figure 2c illustrates the so-called tracking pulses 
where E p and .Fboi track each other over the whole pulse, albeit with a slight shift in time 
(Ford et al. 1995). 

One essential difference in the BR01 study compared to the earlier ones is the use of 
a bolometric flux measure (see also Ryde, Borgonovo, & Svensson (2000)) in which the 
vF v value at the peak is used as a measure for this flux. This method was shown to be 
better than integrating the spectral flux over the BATSE band, as long as the peak of the 
vF v spectrum is in the BATSE window and the power-law indexes, a and f3, do not vary 
significantly throughout the pulse. BR01 also showed that the average distribution of the 
HIC index is the same whether the flux was integrated over the observed band or estimated 
by the peak of the vF v flux, indicating that the bolometric correction for the energy flux in 
most cases does not alter the outcome. The necessity of including a bolometric correction is 
more important for the photon number flux than for the energy flux and is especially true for 
spectra with soft low-energy photon distributions. The usefulness of this method is further 
developed in Borgonovo et al. (2002). The observational results and data used in this paper 
rely on the BR01 approach, so that the discussion will be instrument-independent as long 
as one can find a bolometric correction for the data from the instrument used. 
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Fig. I. — Examples of GRBs (squares and histograms) which follow the relations in equation 
(1) and (2) (dashed lines), (a) the HIC relations and (b) the pulse shapes for the decay phases 
of GRB 920525 (trigger 1625) with rj = 2.4 ± 0.5 and d = 2.0 ± 0.4 and of GRB 950818 
(trigger 3765) with rj = 2.2 ± 0.3 and d = 2.1 ± 0.4. 



-7- 




Fig. 2. — Hardness-intensity correlations (HICs) for four BATSE triggers with 77 indexes 
deviating significantly from r] = 2.0. (a) GRB 921123 (trigger 2067) r] = 2.60 ± 0.15 (b) 
GRB 941026 (trigger 3257) r] = 2.80 ± 0.20 (c) GRB 970420 (trigger 6198) r] = 2.38 ± 0.15 
(d) GRB 980306 (trigger 6630) rj = 1.49 ± 0.04. (see BR01 for details). The long-dashed 
lines show these power-law fits to the pulse decay-phase HICs. Note that the rise-phase HIC 
is included in the figures. The short-dashed lines show the rj = 2 power-law showing that 
the later phases of the pulse in (a,b,c) obey this relation (see discussion in §5). 
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2.2. Pulse Decay Shape 



Investigations of the light curve (over the whole, observed spectral range) during pulse 
decay phases by Ryde & Svensson (2000) have shown that the observed photon flux, F(E)/E 
or F(u) / v, can often be described by a reciprocal function in time. Ryde, Kocevski & Liang 
(2002) derive the corresponding decay shape for the bolometric energy flux (F bol ), which in 
general can be described by 



There is, however, an ambiguity in fitting such a light curve to the data. The values of the 
power-law index d and the time constant r are coupled and, consequently are often not well 
constrained by the fitting. Furthermore, a subjective judgement must be made in choosing 
the transition moment between the rise and the decay phases. As a result the pulse shape 
fitting does not provide as clear a signal as the HIC relation. This ambiguity is illustrated 
by F. Ryde, D. Kocevski, et al. (in preparation), who revisit the BR01 sample and analyze 
the deconvolved light curves and fit the light curves with equation (2). The parameters 
are indeed unconstrained in approximately half of the cases. Two of the constrained cases 
are shown in Figure lb. Kocevski & Liang (2001) and Ryde, Kocevski & Liang (2002) 
introduce an approach to overcome this ambiguity by defining various analytical shapes for 
the whole pulse that include the rise phase (whenever present) and asymptotically approaches 
equation (2) in the decay phase. Kocevski & Liang (2001) studied a sample of 22 pulses 
with good rise phases and found d = 2.2 ± 0.7. 

A corollary of the above two relations is that E p also decays following the above form 
with a different exponent. 



In this paper we will mainly focus on the observed behavior of the spectral evolution as 
described by the HIC and the shape of the pulse decay phase. There are, however, several 
other approaches that can be used to present the observed evolution. One way is to study 
the change in pulse shape and width over different spectral bands and the time lag between 
these bands (see, e.g., Norris, Marani, & Bonnell (2000)). To illustrate this we present in 




(2) 




(3) 



2.3. Change of Pulse Shape with Spectral Band 
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Figure 3 the light curves of the four cases shown in Figure 2. For each case we show the 
light curve in two energy bands. For more detailed discussion on this approach we refer to 
F.Ryde, D.Kocevski, et al. (in preparation). 



3. TIME SCALES 

There are many time scales in a GRB. In this paper we are interested in the effects of 
the curvature of the emission front and the time scale it introduces in the pulse shape. In 
order to find the importance of this effect its time scale should be compared to the radiative 
cooling time and the dynamic time for the crossing (or merger) of two shells. 

Estimation of the cooling time, which we assume is comparable or longer than the 
acceleration time, requires a knowledge of the emission and energy loss processes of electrons. 
The details of the acceleration mechanism are not known, but it is generally assumed that 
the dissipation of the kinetic energy of the fireball gives rise to a power law distribution 
of electrons. The minimum random Lorentz factor of the electrons is proportional to the 
relative Lorentz factor; 7 e = £ e (m p /m e )(r rc i — 1), where m e and m p are the electron and the 
proton masses, respectively, and r rc i = ^^(l — P1P2) ~ Fi/(2r 2 ), which is of the order 
of 10 or lower, is the relative Lorentz factor between the two interacting shells (r x > T 2 ). 
£ e , and £ B introduced below, are the fractions of the post-shock, random energy density 
that resides in the electrons and the magnetic field. The magnetic field strength is given by 
B = (327r^B^ / ^pC 2 r rc i) 1 / 2 ~ 2 x 10 4 Gauss, where we assume equipartition and a comoving 
electron density n' = Lk in /47i\R 2 r 2 m p c 3 ~2x 10 10 cm~ 3 , for an assumed kinetic luminosity 
of L kin ~ 10 53 erg/s, distance R of 10 15 cm, and a shock compression ratio of 4. 

In the rest frame of the progenitor the cooling time scale for a particle of comoving 
energy E' and loss rate E' is dilated to ri oss = TE'/E'. As the emitting material is moving 
towards the observer the rest-frame time scales will be compressed by a factor (1— /3) _1 ~ 2r 2 
in the observer frame (see below). If the pulse shape is determined solely by the radiative 
cooling process then, for the above parameters and assumptions, the observed decay time 
scale of the pulse will be 



E 2>Tim e c 
Tc °° l = 2T\E S + E 1C \ = rV^B^I + U iad /U B ) 

~ 6 x 10- 5 s(r rcl - I)- 2 ( V— V% V 1 (4) 

v ; V 10 5cm / V 100 / V 10 / 

Here subscripts s and IC refer to the cooling rates due to synchrotron and inverse Compton 
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losses in radiation and magnetic field energy densities of C/ ra( j and Ub = B 2 /(8n). This is a 
much shorter time scale compared to the dynamic and curvature times discussed below. 

The curvature time scale arises from relativistic effects in a sphere expanding with a 
high bulk Lorentz factor T. Due to the curvature of the shell there will be a time delay 
between the photons emitted simultaneously in the comoving frame from different points on 
the surface. Figure 4 shows the geometry of the situation. Due to the relativistic aberration 
of light, isotropically emitted radiation in the comoving frame will be beamed into a cone 
with opening angle 9 ~ Only photons emitted from the fireball surface within a narrow 
cone of opening angle ~ T^ 1 around the LOS will be detectable by the observer. The typical 
time delay is thus r ang = R[l — cos(l/T)]/c which for large T is approximately R/(2T 2 )c. 
This gives a lower bound for the observed duration of a pulse: 



T - = L7s (lo^)(^)" 2 ' (5) 

The dynamic time scale for a single pulse is the actual crossing (or merger) time of one 
shell with another. Often the shell collision is assumed to be an inelastic collision and the 
merged shells expand as a single shell. The shell crossing time is r^ yn = A'/v' sh , where v' sh is 
the velocity of the shock in the comoving frame of the preshocked flow. The initial value, and 
the evolution with radius, of the shell width A' are not well understood and depend, among 
other factors, on the structure and internal dynamics of the shell and on its interaction with 
the external medium. If the Lorentz factor is constant over the shell, or if the shell is confined 
by some mechanism, then the shell width will be independent of radius and time, A = A . 
But if there exists a differential flow with a faster leading and a slower trailing edge with 
velocities f3 T and f3 B , and corresponding Lorentz factors T r and T s , respectively, then 



A(R) = R((3 T -(3 S ) 



R 

s 



1 



R 



where for the last relation we have assumed that r s C T 
R/(2T B ). This relation will be applicable at radii R > R^ = 
exceed the initial width. The shell crossing time is then r^ yn 



2r 2> 

s 

In the comoving frame A'(R) 



(6) 



A r 2 , where the spread will 
= A'/v' h = R/(T s v' sh ). In the 
observer frame these times will be time-dilated and affected by the motion of the fireball 
towards the observer: 



r dyn 



R 



'dyn 



2r 4r s ry 



sh 



sh 



R 



100 



-2 



(7) 



For the assumptions described above the dynamic time scale is somewhat shorter than (and 
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has a similar dependence on R and T) as the curvature time scale; r ang ~ 2(r s /T)rd yn . 
However, at smaller values of R < Rb, where the spreading is negligible, the dynamic time 
scale may exceed the angular one by Rb/ R. Since, in general, the pulse width is proportional 
to R, the latter situation is more likely to arise in shorter pulses. 

Figure 5 shows the radial dependence of the three time scales, r coo i(eq. 4), r ang (eq. 5), 
Tdyn (eq. 7), for two values of the bulk Lorentz factor T, for T ~ T s , and for Rb = 10 11 
cm. From this we conclude that in situations with sufficiently narrow shells, the curvature 
effects can indeed be the determining factor for the pulse decay time scale. Correspondingly, 
if these effects are shown to be important in the observed pulses, this could set a constraint 
on the thickness (value of r s /r r ) and the spreading of the shell (values of Rb and r s ) and/or 
on the distance R where the radiation is emitted. The curvature effect becomes important 
if Tang > T dyn . This happens for shell widths A' ^ [3' sh R/T. Even in the case with linearly 
increasing shell-width the angular effect will be noticeable and can be totally dominating if 
the shell is narrower (e.g. if the shell- width stays constant). Spada, Panaitescu, & Meszaros 
(2000) simulated internal shocks and found that with a linear shell broadening, as above, 
the angular spreading and the shock crossing times are comparable. 

In the following we consider a scenario where the rise phase of a pulse corresponds to the 
merging phase or the dynamic evolution, while the decay phase arises as a purely kinematical 
effect due to the curvature of the relativistically expanding shell. The acceleration and 
radiative cooling times are assumed to be much shorter than or at most comparable to the 
dynamic time. 

4. RELATIVISTIC KINEMATICS 

The radiation from a relativistically expanding plasma sphere (fireball) will have a 
unique signature in the observer frame. We start by studying the radiation observed from 
an infinitely thin spherical shell and thereafter generalize the problem to broader shells. We 
assume that the visible part of the fireball is spherically symmetric and homogeneous, i.e., 
only one parameter is enough to characterize the properties of a patch of emission. This 
type of problem has previously been discussed by, among others, Fenimore et al. (1996), 
Granot et al. (1999) and Eriksen & Gr0n (2000). 
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GRB 921123 GRB 941026 




seconds since trigger 



Fig. 3. — Background-subtracted, count light-curves (64 ms resolution) for the four triggers 
in Figure 2. For each trigger the light curve for two different spectral bands are shown; in 
gray, BATSE channel 2 (55 - 110 keV) and, in black, BATSE channel 3 (110 - 325 keV). The 
time interval used for the HICs in Figure 2 are indicated by the dashed lines. 




Fig. 4. — Geometrical sketch of the visible part of the fireball. Due to the relativistic 
aberration of light, the main part of the radiation comes out within a small cone with half 
opening angle 9 ~ Note that the line of sight (or LOS) and the jet axis do not need to 
be parallel. The photons from A are boosted by 2r but those from B, which are delayed by 
At = R/c(l — /i), are boosted by [r(l — where \x = cos9. 
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Fig. 5. — The radial dependence of the various time scales. The solid lines are for a bulk 
Lorentz factor of T = 100, while the dashed lines are for T = 1000. For the cooling time 
we have used eq. (4) with r rcl = 10 and a kinetic luminosity L kin = 10 53 erg/s, and the 
angular time is given by equation (5). For the dynamic time we have used equation (7) with 
T s ~ T = 100 and (3' sh = 1. If the shell width is constant, A = Ao, then there will be a 
minimum dynamic time scale at all R and be represented by a horizontal line in the figure; 
Assuming A = 10 9 cm, r dyn = 1.6 x 10" 4 s ^ h _1 (r/100)- 1 (A /10 9 cm). 
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4.1. Lorentz Boosting and Spectrum Evolution 

The Lorentz boosting factor for transformation from the comoving frame to the observer 
frame of photons emitted into the LOS from different locations on the surface, defined by 
the angle = arccos \x shown in Figure 4 is 

For small angles this reduces to T>(T) = 2T, which is the boost factor used in discussing 
GRB emission, when the angular dependence is not considered important (e.g. flat shell 
perpendicular to LOS). 

If we set = at the point where the flow velocity is parallel to LOS (point A in Fig. 
4), then the difference in light travel time At between photons emitted along the LOS from 
this point and a point at an angle 6 (point B in Fig. 4) is At = Rq{1 — fx)/c , which gives 
fx — 1 — cAt/R . Inserting this into equation (8) we find 

p(r - At) = r ( iVftW (9) 

For highly relativistic outflows (3=1 — (2r 2 ) -1 , and 



V(T, At) = 2T(1 + At/T ang )-\ r ang = R /(2T 2 c). (10) 

An obvious outcome of this is that if the emitted spectra from different parts of the shell 
are identical, then the observed spectrum will be boosted (gradually redshifted) in time by 
a factor V as different parts of the surface come into view. In particular the peak energy 
will evolve as 



where E' is the peak energy in the comoving frame (which is the same at all angles), E = 
E' p 2T, and we have set At = t, which assumes that t = for radiation observed from point 
A in Figure 4. 



4.2. Energy Flux and the Light Curve 



We want to determine the pulse shape, in bolometric energy flux, assuming it to be 
entirely caused by the above angular dependence of the Lorentz boost variation. This means 
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we are considering a situation where Td yn and r coo i <C r ang . We, therefore, initially assume a 
simple isotropic volume emissivity. 

f = Z'5(t' -t )5(R' -R ), (12) 

where £' is the total, rest frame, energy surface-brightness. More complex temporal and 
spatial distributions will be discussed in the next section. [We follow the formalism outlined 
in Granot et al. (1999) and refer to Figure 4.] We also introduce the following variables: 
the luminosity distance d L , the viewing angle a = sm6R/d L = (1 — fi 2 y^ 2 R/d L so that 
ada = (R/di,) 2 fidfi. Using the Lorentz invariance of j/u 3 (Rybicki & Lightman 1979) we 
can write the observed intensity (neglecting cosmological redshift factors) as 



J jds = J j'V 3 dR/fi (13) 



and the observed bolometric energy flux as 



bol 



J fildtt = J J J j'V 3 dRsmadad® (14) 

where ds is the element of length through the shell along the LOS, $ is the azimuthal 
angle and dVL is the solid angle of the source as seen by the observer. Assuming azimuthal 
symmetry around the jet axis ( f d$ = 2n) we get 

F ho i(t) = 2J2 R2dR / vdfxV 3 (4nf). (15) 

L J J — 1 

Using the transformations \dR/dR'\ = \dt b s /dt'\ = X> -1 , and that dt = (i?/c)rf/x 
we get S(t' - t )5(R' - R ) = c (RTV)' 1 6(fi - pL )S(R - R ). The observed flux then can be 
written as 

F ho i(t) = F 0f iV 2 (fi)/T 2 , with F = 2nZ / R cT/dl (16) 

where the time dependence arises from the relation (j, — 1 — ct/R . For a highly relativistic 
outflow (where 9 max ~ 1/r) the projection factor [i — cos 9 in equation (16) is less than 
1 — (2F 2 )" 1 and can be ignored so that F ho \(t) = F V 2 (fi)/T 2 . Note that the above result is 
independent of whether the outflow is completely spherical or confined or collimated into a 
jet of opening angle 9j et , as long as 9- ]Ct > T^ 1 . 

Now with the help of equation (9), and identifying the time delay At with the time t as 

above, we get 

i 71 boi(t) = F /(l + t/r ang ) 2 . (17) 
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This surprisingly is very similar to the observed behavior of the small sample of well-isolated 
pulses described in §2. 



4.3. Flux-Spectrum Correlation; HIC 

In the 5-function approximation, and if the bulk Lorentz factor and if the rest frame 
flux and spectrum are independent of position along the surface of the shell (i.e. we have a 
homogeneous source), then equations (11 and 17) will describe the spectral and flux evolu- 
tions adequately. This can be translated into a measurable hardness-intensity correlation if 
we represent the hardness by the value of E p . The observed HIC (eq. [1]) follows directly: 
F oc E£, with a power- law index of rj = 2, which is equal to the observed average value found 
in BR01. 

Another way to look at this relation is through the time evolution of the energy fluence 
defined as E(t) = J* F(t)dt. Using equation (17) it can easily be shown that 

S(t) = F r[l-(l + t/r ang )- 1 }, (18) 
which when inserted into equation (11) gives 

E p = E pfi -^, (19) 

which is the so-called hardness- fluence correlation (HFC). The inverse of its proportionality 
constant is usually denoted by <3> (Liang & Kargatis 1996; Ryde & Svensson 2000), so in 
this description $ = F T ang /E Pt0 . 



5. FINITE DYNAMICAL TIME 

The above analysis assumes a very short dynamic, as well as cooling, time scales relative 
to the curvature time delay. This allows the delta function time profile approximation used 
in §4.2. However, as mentioned above and shown in Figure 5, the dynamic time rd yn , which 
in our model defines the rise time of the pulse, could be comparable to (and in short pulses 
it may exceed) the curvature time scale. Therefore we must examine the effects of finite r dyn . 
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5.1. Bolometric Light Curves 

For a shell shining continuously for a finite period of time the resulting pulse profile 
will be a convolution of the comoving frame emissivity time profile, say j'(t') = H'f(t'), 
with f f(t')dt' = 1, and the profile due to the curvature effect (eq. [11]) obtained from an 
impulsive intrinsic emission profile, / = S(t' — t' ). The emissivity profile, when transformed 
to the observer frame, retains its form but its characteristic time scale is scaled by the boost 
factor; r dyn = r dyn V. The resultant pulse profile then becomes 

/+oo 
f(t-x)dx/(l + x/r ang ) 2 (20) 

Figure 6 illustrates the effects of the finite emission time scale for the intrinsic pulses. 
In the three top panels the intrinsic pulses were modeled as decaying exponentials with time 
scale r dyn . The shape of the convolved light curve is determined by the ratio of the intrinsic 
emission time scale r ang to that of the curvature Td yn ; 7£ = r ang /rd yn . In general the convolved 
pulse will resemble a FRED (fast rise and exponential decay). Different ratios of the time 
scales give rise to pulse shapes that encompass the actual observed pulse shapes by BATSE. 
For long duration pulses the convolved light curves will asymptotically reach the d = 2 
form of equation (2) and manifest the curvature effect. But for short and weak bursts this 
stage may not be observable and one sees a pulse shape determined by the dynamics of the 
shell crossing. The general behavior is well illustrated by these three cases. However, other 
intrinsic pulse shapes do give rise to some variations. For instance, including an intrinsic rise 
phase will broaden the rise phase of the convolved light curve, as illustrated by the bottom 
panels in which the intrinsic pulses were modeled by symmetric Gaussians of half width r dyn . 

It should be emphasized again that in the discussion above, and in equation (20) in 
particular, we are dealing with the bolometric flux, using the approach introduced in BR01, 
who represent this flux by EF E evaluated at E p . In general, however, for a finite observation 
band there will be a time dependence term beyond what is shown in equation (20). If, over a 
finite dynamic time, the spectrum evolves, i. e. the E p changes, the flux contribution within 
the observed band, and consequently the bolometric correction will change (the amount of 
changes depending on the values and/or evolution of indexes a and 0). This will affect the 
above light curve as well as the HIC relation described next. In an approach different from 
BR01, to account for this issue, Fenimore & Sumner (1997); Fenimore et al. (1996) instead 
assumed a photon number spectrum described by a single power law with index equal to the 
averaged value (a + f3)/2, and calculated the observed photon flux over a certain band pass 
which thus changes in the observer frame with time. The relation between the photon flux 
and the peak energy found by this method was studied by Soderberg & Fenimore (2001) 
who concluded that a long cooling time was necessary to explain the data. 



Time (seconds) 

Fig. 6. — Light curve behaviors in the observer frame, with varying parameter combinations 
of the exponential decay constant Td yn and the angular time scale r ang . The dotted, dashed, 
and solid lines represent the comoving light curve, the relativistic angular spread function and 
the convolved light curve. Top panels for a one-sided (zero rise time) exponential intrinsic 
light curve / oc exp(— t/r dyn ), t > 0. Note that the convolved light curve has a finite rise 
time. The lower panels for a symmetric Gaussian intrinsic light curve, / oc exp(— t/rd yn ) 2 , for 
— oo < t < oo. The comoving light curves and the angular spread functions are normalized 
to have the t — at the peak of the observed (convolved) pulse. 

5.2. Spectra and HICs 

The effects of a finite emission pulse on the spectral behavior are more complicated. 
This is because the instantaneous spectra, as seen by the observer, will be a superposition of 
spectra from several annuli emitted at different times having different fi = cos 6, boost factor 
T>(fi), and perhaps different spectral hardness (say different E p ) and shape (different indexes 
a and j3). Let us assume that the intrinsic spectrum, as observed in the rest frame, can be 
described by two spectral indexes a and (3 and a break or peak energy E p , say a Band et al. 
(1993) function, B a ^(E/E p ), with J °° B dE = 1. To simplify matters, we further assume 
that a and j3 are constant so that the only spectral variation is due to changes in E p . We 
will represent the evolution of E p by an intrinsic HIC relation with index 77 int , F oc E p mnt , 
which means that E p (t) = E Pj0 [f(t)] < -~ 1 / r ' int \ The observed flux spectrum at any time can 
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thus be written as 



F(E,t)=F 



oo 



(1 + x/r ang ) 



fit - x) 




( 



E p (t,x) 



E 



dx, 



(21) 



with 



E p (t, x) = £ p ,o(l + x/r^-'Uit - x)}^ 1 ^ 



(22) 



When integrated over all energy, equation (21) reduces to equation (20). Note that the 
resultant spectra or light curves depend primarily on the ratio of the time scales, TZ = 
Tang/Vdyn defined above. Secondary factors are the intrinsic pulse shape and HIC relation. 

Before showing the resultant HIC for different cases we first discuss the details of 
the formation of the observed spectra. For the purpose of illustration we assume a = 
— 2/3, /? = —2.5 and three different values for the index 77 int . To simplify the description 
of how the contributing spectra make up the integrated spectrum we start with the case 
of an intrinsic light curve with an abrupt rise phase and a simple exponential decay phase: 
f(t) = Tdyn -1 exp(— t/r dyn ), for t > 0. The effects of intrinsic pulses with a finite rise phase 
are discussed below. Note that an intrinsic pulse with only a decay phase will still produce 
a rise phase for the observer due to the convolution described above. For such an intrinsic 
pulse, the upper limits of the integrals in eq. (20) and (21) are zero, so that at any given 
time t the observer receives signals from annuli extending from 6 = to 6 = 9 max (t) or 
1 < pL < /i min = (1 — fict/Ro). Different angles 9 sample different stages of the intrinsic 
spectral evolution due to the angular spreading of the signal, with the spectrum from # max (t) 
being that emitted at the beginning of the pulse and the spectrum from 8 = reflecting that 
of a later time in the pulse. Equation (21) gives the convolved spectrum as a superposition 
of different spectra from different angles or different times. 

In the absence of the curvature effect (i.e. for a flat shell, Rq — > oo, or r ang and TZ — > 0) 
one would observe the intrinsic HIC relation. It follows then that for a finite r ang but for 
TZ <^ 1, this is the HIC relation that will be observed except for a short time t ~ r ang at 
the beginning (and the rising phase) of the observed pulse. In the opposite limit of TZ ^> 1, 
the delta function approximation result with the HIC relation F bo i oc E p is obtained, again 
except for a short time t ~ r dyn at the beginning (and the rising phase) of the observed pulse. 

For TZ of the order of unity, one expects a combination of these two behaviors. For the 
assumed exponential intrinsic light curve, which (for TZ ~ 1) dies more quickly than the light 
curve shape induced by the curvature effect, we expect an initial phase when the HIC obeys 
the intrinsic form (rj ~ ?7 int ) followed eventually by the relation expected from the curvature 
effect (rj — 2). This transition can be seen in Figure 7 for the three cases with TZ = 0.5, 1 and 
2, from the top to the bottom panels. In each panel we show total spectra for four different 
observed times t (solid curves). The circles show the peak flux and the photon energy at the 
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peak (i.e. the observed E p ) for these four and several more total spectra. This is what the 
observed HIC will look like. The light dotted lines show the assumed intrinsic, ?7 irit = 1, and 
the curvature-induced, r] — 2, HIC relations. Each one of the total spectra is a sum of the 
spectra from different angles < 9 < 9 max (t) with different E p s having suffered different 
boosts. The peak flux and the peak energy of spectra from several angles are shown by the 
crosses for each of the four observer times. The crosses furthest to the left represent the 
9 = spectrum which is identical to the intrinsic spectrum (of course boosted by 2T 2 ) and, 
therefore, lies on the dotted line with slope i] int = 1. As we move to the right we see spectra 
from successively larger 9s (with varying boosts) till the crosses furthest to the right which 
are for # max , i.e., the radiation emitted at t — 0, with the flux and E p boosted according the 
relations described in the previous section, and lies on the dotted line with the slope r\ = 2. 

For the lower panel, TZ — 2, the most intense radiation comes from # max . Consequently, 
the total spectra obey the large TZ, or delta function, HIC relation (rj = 2), except for a 
relatively short but noticeable initial phase. For smaller TZ this initial phase increases in 
length. As evident from the relative position of the crosses in the panel for TZ = 0.5, for the 
first two times shown, the emission is dominated by the radiation from 9 = and we obtain 
a well defined HIC with rj = i] int = 1. But by the time corresponding to the fourth total 
spectrum the emission from 9 max dominates and the HIC reverts to the r\ = 2 case. The 
location of the crosses and the final HIC relation for the TZ = 1 case is intermediate between 
the above two cases. 

In order to further demonstrate the relative effects of the intrinsic and curvature-induced 
HIC relations, in Figure 8 we show two cases with TZ — 1 but with intrinsic power law indexes, 
77i nt = 0.5 and rji nt = oo. The latter corresponds to a constant E p independent of the flux. In 
Figure 9 we show only the resultant HIC relation for five different values of TZ for the above 
three intrinsic HIC indexes r] int = .5, 1 and oo. 

Several conclusions can be drawn from these figures. The first is that the observed 
evolution of E p may have nothing to do with the intrinsic evolution of the emission process, 
so that care is necessary in the interpretation of the observed spectral parameters. This is 
most clearly evident in Figure 8b where a non-evolving spectrum appears to have a HIC 
index of rj = 2. It is also clear (from all of the figures) that this most commonly observed 
HIC index does not require a very large value of the ratio TZ. It is established after a few 
dynamic times for TZ = 1 and is omnipresent for TZ > 2, independent of the value of r] m f On 
the other hand, the intrinsic HIC is essentially what is observed (at least when the flux is 
high, say down to 1 % of peak flux) for TZ < 0.5. 

It is clear that if there is some dispersion in the intrinsic power-law index, i]i nt , for 
instance, if it has a broad, uniform distribution between and oo, then the observed dis- 
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Fig. 7. — Observer frame spectra (solid lines) at four different times (1, 2, 6 and lis; where 
the first light is detected at t — s). These spectra are found by integrating the intrinsic 
spectra over 9. The positions in the F — E p plane of the peaks of some of these contributing 
spectra are indicated by large crosses. The upper and lower dashed lines represent, respec- 
tively the locii at = which follow the intrinsic evolution of the HIC relation (here assumed 
to be a power law with intrinsic HIC index, rj- mt = 1.0), and at 6 — max obeying the HIC 
relation expected from the curvature effect (with a power-law with index r\ = 2). The circles 
show the HIC relation of the convolved spectra in time steps of 1 s, up to t — 60 s: (a) 
ft = 0.5 (b) Tl = 1.0 (c) K = 2.0. 
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Fig. 8. — Same as Fig. 7 for 1Z — 1.0, but two different intrinsic HIC relations (a) r]i nt = 0.5 
and (b) f]i nt = oo corresponding to a constant E' At). The softening of the observer frame 
low-energy spectra is clearly seen in (a) while (b) demonstrates that for steep intrinsic HICs 
(?7i n t > 2), the slope of the observed HIC will mainly be 2, independent of the intrinsic slope 
(see also Fig. 9a). 
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Fig. 9. — Observer frame HIC relations, for different values of TZ, produced by three different 
intrinsic HIC relations: (a) i]i nt = oo (b) i]i nt = 1 (c) i]i nt = 0.5. In panel (a) the time steps 
are 2 s and 1Z = 0.05 (crosses), 0.1 (x-es), 0.2 (diamonds), 0.5 (inverted triangles), f.O 
(triangles), while in panel (b) and (c) the time steps are 1 s and 1Z = 0.05 (crosses), 0.5 
(x-es), 1.0 (diamonds), 2.0 (inverted triangles), 20 (triangles). Note the presence of various 
concave HICs which can explain observed bursts with index > 2, when fitted with a single 
power law. 



-24- 



tribution will be peaked around r] = 2. If the intrinsic source of radiation is synchrotron 
emission the preferred value of i] int will lie between 1 and 2 (Lloyd & Petrosian 2002) and 
the observed distribution is expected to be even more peaked around 2. However, a variety 
of other possible values of this index and deviations from a simple power law will also be 
present. This is true especially for values of 1Z ~ 1. In addition, limitations imposed by 
the instrumentation make it possible to follow only a portion of the HIC. This is because 
the signal-to-noise of the observations limits the dynamic range of the observable flux and 
the finite band width of the instrument introduces bias against spectra with E p outside the 
band. As is evident from Figures 7, 8, 9 some of the curves do not follow a perfect power 
law but are concave, which is even more pronounced for smaller i] int . Thus a simple power 
law fit to these curves, excluding the few early phase points which may appear as part of 
the rise phase of the pulse, will yield a steep HIC with rj > 2. In fact inspection of the 
pulses with large measured rj- values in the BR01 sample and illustrated in Figure 2 show 
some resemblance to these concave curves. 

The second conclusion is that because the observed spectrum is a superposition of many 
intrinsic spectra, it will necessarily be broader than the intrinsic spectrum and its spectral 
parameters could have complex relations with the intrinsic ones. Two obvious effects of this 
broadening can be seen in the above figures. One of these is that the value of E p is less 
well determined because of the flat tops. This is the cause of the sharp transition from a 
high (intrinsic) to a low (boosted) value of E p for the cases with 1Z < 1, and 77 int = oo. The 
second effect of this is that for fits to a finite spectral range of observations the slope of the 
flat portion will be identified mostly with the low energy index and result in a = — 1 instead 
of assumed a = —2/3. In some cases this portion could be identified with high energy index 
and yield (3 = — 1 instead of —2.5. This effect is more pronounced for smaller values of ?7i n t 
(compare Fig. 8a and Fig. 7b) because E p undergoes a larger variation for small changes 
in the flux. More than half of all pulses observed by BATSE have spectra that do change 
their shapes in time, often showing a low energy softening (Crider et al. 1997). The above 
figures show that such behavior can be produced purely by the curvature effect without the 
intrinsic spectrum changing. 



5.3. Effects of Different Intrinsic Light Curves and Rise Times 

In the discussion above we have used a simple prescription of the intrinsic light curve 
with only a decay phase: exp(— t/rd yn ). Its actual shape is revealed to the observer only 
if 1Z > 0.5. Otherwise any claim about the pulse shape is pure speculation. We have also 
explored other shapes and found that the resulting HICs show the same qualitative behavior, 
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albeit that a given behavior is found for a different value of the characteristic time scale, 
Td yn . These similarities and differences can be seen by comparing the relations marked by 
squares in Fig. 10a, derived for a pulse with a Gaussian decay phase [fit) oc exp[— (t/r dyn ) 2 ], 
t > 0, solid line], and for an exponential decay [dashed line; compare Fig. 8a]. 

A second aspect of the intrinsic light curve which must be explored is its finite rise time 
when this time is not much shorter than T& yIl and/or r ang . The curvature effect during the 
rise phase of the intrinsic pulse is somewhat different from its effect during the decay phase 
described above. Figure 10b shows this for an intrinsic, complete, not one-sided, Gaussian 
light curve. During the intrinsic rise phase (points 1 and 2 in the figure) the flux will always 
be dominated by the emission from 9 = 0. This is because (i) larger angles sample earlier 
times in the intrinsic light curve (which are weaker) and (ii) the boost factor becomes smaller 
at larger angles. Therefore, at the early time of the observed rise phase (before the intrinsic 
decay has started) the H1C will simply follow the intrinsic 77^, largely independent of 1Z. 
During later times, when the intrinsic decay starts to be visible at the smallest angles, the 
behavior outlined in the previous section will occur. The rise phase will only affect the 
observed HIC in the early part of the rise phase. This is also shown in figure 10a where the 
HICs are plotted for a one-sided Gaussian with only a decay phase (squares, solid line) and 
and a symmetric Gaussian (crosses). Including the rise phase will produce an up and down 
(soft-hard-soft) or so-called tracking pulse while the pure decay will produce a hard-to-soft 
evolution, discussed above and by Ford et al. (1995). 

In summary, for the intrinsic rise phase the flux will be dominated by 9 = and the 
HIC will follow the intrinsic i]- mt . This makes the rise phase of a tracking pulse of interest as 
it reveals the intrinsic HIC. For the decay phase the ratio 1Z determines the behavior. For 
1Z >> 1 the flux from 6> max dominates and rj — 2, while for 1Z <C 1 the flux is dominated 
by 9 = and the HIC follows r] int . For cases in between, a HIC with 77 = 2 dominates, for 
most of the observable part of the pulse, down to 1Z ~ 0.5, the exact value depending on 
the intrinsic light curve shape. In general the decay phase HIC follows initially a concave 
curve sometimes resembling an S-shaped curve. The transition in these cases is similar but 
much broader than those referred to as track jumps by BR01 (the case with two parallel HIC 
curves displaced from each other), which was explained as two separate spectral components 
coming from two overlapping pulses. The S-shape feature, with a transition from i] = r] int to 
i] = 2, can explain the more moderate changes seen in Figure 2, and demonstrates that the 
HIC shape can be used as a diagnostic for the value of 1Z and the intrinsic light curve. 
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Fig. 10. — Variations in the spectral behavior due to the variation in the intrinsic light curve, 
(a) Same as figure 9 but for a Gaussian intrinsic light curve and 1Z = 0.25, a one-sided, pure 
decay Gaussian produces a so-called hard-to-soft pulse (squares), while a symmetric Gaussian 
(crosses) produces a so-called tracking (hard-soft-hard) pulse. For the latter case the rise 
phase will follow the intrinsic HIC. For comparison, the dashed line shows the HIC for a 
one-sided exponential decay with 1Z = 1.0 from Fig. 9b. (b) Same as Fig. 7 but for a 
Gaussian intrinsic light curve with equal rise and decay time scales. Circles show the HIC 
relation and the large crosses show the peaks of the contributing spectra. Times 1 and 2 are 
during the rise phase while 3 and 4 are for the decay phase. 
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6. MODIFICATIONS OF THE BASIC MODEL 



In deriving the above equations and the HIC we have made several simplifying assump- 
tions. We therefore briefly discuss how the results are affected when we relax some of these 
restrictions. 



A limitation of the results presented in §3 is the delta function representation of the 
spatial distribution of the emission. The actual emitting region most likely has a finite 
width, less or equal to the width A' of the shell. Here we assume it to be equal to the 
shell width, with f oc g(R',A') instead of the delta function description used in §4. The 
above approximation is valid for A'/R <C 1. This can be demonstrated assuming a Gaussian 
distribution, g(R',A') = 2/(A'y/n) exp [-((#' - R)/A') 2 }. If we replace the 5(R'-R) factor 
in equation (12) by this form and carry out the derivation again, then equation (16) becomes 



If the relation between A' and R is the one given by equation (6) this becomes F(t) ~ 
F fi(V/T) 2 (l + 1/r), which shows that the width of the spatial distribution affects the flux 
(negligible for T 1) but has no effect at all on the spectrum of the emission. 



Another tacit assumption has been that the emission surface is spherical, either forming 
a complete sphere or, if it is confined into a jet of opening angle 9 Jct) that the LOS intersects 
the jet at an angle more than 1/r from the edge of the jet. If the LOS is close to the edge of 
the jet then an extra steepening of the HIC can arise, apart from the nominal value for the 
angular spreading case (77 = 2). Let us consider this effect for the delta function emission, 
or 1Z — > 00. During the early stages of the pulse an observer will see full (2tc) annuli with 
E p and F varying as described in §4. However, when the azimuthal symmetry is broken, 
i.e. the annuli are no longer whole, the flux integrated over the partial annuli will be less 
compared to what is expected from a full annulus. However, the observed value of the E p 
will not be effected as the spectrum will not change. This will lead to a break in the HIC, 
from i] = 2 to a larger 77. A corresponding behavior will also be seen in the light curve. This 
could explain the most extreme cases in the BR01 sample. 



6.1. Spatially Broad Emitting Region 
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6.2. Geometrical Effects 
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6.3. Inhomogeneities 

In the description of the relativistic outflow we have assumed that the angular distribu- 
tion of the Lorentz factor is homogeneous, i.e., it is constant over the angles and and that 
at angles larger than 6- ]et) Y drops abruptly to zero. However, this is just a simplistic descrip- 
tion and does not necessarily describe the actual situation. Deviations are most probably 
important at the edges of the jet. Therefore, if the beaming angle is of the same size as the 
jet opening angle and/or if the LOS is close to its edge, then the intensity will be affected. A 
variable Lorentz factor will give rise to a slight deviation from spherical symmetric outflow. 
This will affect only the shape of the light curve and the time evolution of the observed E p , 
but not the HIC because both F and E p will be affected same way. 

6.4. Spectral Variations 

We have also not considered the effects of the intrinsic, rest frame, spectral-shape varia- 
tion (changes in a and (3) during the pulse, but have rather concentrated on the exploration 
of the variation of E p and flux. As shown above, there is a softening that occurs naturally 
in this scenario and hence could be one of the reasons for the observed softening without 
the intrinsic spectrum necessarily changing. For a more general study, individual cases need 
to be examined. This emphasizes again that the interpretation of the observed spectra and 
their relation to the source spectrum, and consequently the emission mechanism, is not 
straightforward. This aspect of the problem is, however, beyond the scope of this paper, but 
will be dealt with in a future paper. 

7. SUMMARY AND DISCUSSION 

We have examined the effect of differences in light travel time due to the curvature of 
the expanding shell and determined to what extent it can affect the width and shape of 
pulses and their spectral time evolution. The energy flux F and the photon energies E will 
be affected by the angle-dependent Lorentz-boost factor £>(//); F oc V 2 and E oc V. The 
peak energy of the vF v spectra, E p , will thus follow a hardness- intensity correlation (HIC) 
F oc E v with i] = 2. Furthermore, the decay phase of a pulse will follow the form (l+t/r ang ) 2 . 
We show that this effect should be important for a reasonable choice of parameters (Lorentz 
factor, burst energy, shell width etc.) and that these characteristics agree with the average 
behaviors found in pulses. 

However, the curvature effect can not alone explain the large observed dispersion of the 
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HIC index r] . Cases with rj largely different from 2 we believe are produced by a finite 
dynamic time, r dyn 7^ 0. The resulting spectral/temporal behavior depends mainly on the 
ratio 7Z = r ang /T^ yn . The intrinsic HIC (assumed to be a power law with an intrinsic index 
Vim) w ih be revealed when 7Z < 0.5 while the behavior expected from the curvature effect 
with a HIC index rj = 2 will dominate for 1Z > 2. For intermediate TZs the F — E p relation will 
deviate from a pure power-law, having a more concave shape (in a log — log plot). A general 
softening of the spectra with time, which has been observed, is also expected, independent of 
any changes in the intrinsic spectrum, and therefore independent of the physical environment 
where the pulses are produced. 

An important conclusion of this work is that one must be very careful in the interpreta- 
tion of the observed light curves and spectra, their parameters and evolution. This is because 
we have shown that the observed light curve will in most cases be different from the intrinsic 
one and the observed spectra will have a complex relation to the intrinsic ones. The spectra 
in the observer frame will be broader and, for instance, the low-energy power-law slopes a 
will be softer than the intrinsic ones. In some cases, flat-topped spectra are produced which, 
in the observer frame, appear to have either a or (3 = —1. Furthermore, we also explain the 
occurrence of pulses whose E p track the flux up, and down, during the rise and decay phase, 
respectively, as well as the occurrence of pulses where the hardness declines monotonically 
independent of the rise and fall of the flux. 

Ultimately, we wish to determine the characteristics of the intrinsic emission, namely 
F'(t), E '(t), i] int , and if possible the spectral power law slopes a and f3. In addition, we 
want to determine the distance R from the progenitor where the fireball emits the 7-rays 
and to discern something about the shell width A' and/or its spreading. Fits to the HIC 
and observed light curve will be able to reveal 7Z and F'(t). Below we discuss the principle 
diagnostics that can be made for three different situations. The value of the bulk Lorentz 
factor T remains an unknown parameter. 

Case I. The observed HIC is a pure i] = 2 power law (e.g. pulses in Fig.l): According 
to our model the curvature effect is dominant (It could, however, also be due to an intrinsic 
rj = 2 HIC). The observed light curve will (asymptotically) follow equation (2) with d = 2, 
from which one can determine the value of r ang . This time constant determines the distance 
at which the shell lights up: 



The observed light curve can, in principle, be deconvolved, with equation (17) as the impulse 
response, to obtain F'(t). A better knowledge of F'(t), then gives a more accurate value 
of 1Z (and thereby Td yn ) from a fit to the HIC. Knowing r dyn one can put constraints on 
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A'/Psh or /3 sh '(T s /T). Furthermore, the observed, instantaneous spectra will be results of 
integrations of the intrinsic spectra along a rj = r] int power law. Combining this knowledge 
with the observations could reveal E' (t) and r] int ~ a obs and possibly a int . 

Case II. The HIC is a pure power law with index 77, substantially different from 2.0 
(e.g. pulse in Fig. 2d): Here, rj = r]- mt and the light curve should reflect the intrinsic F'(t) 
(smoothed somewhat by the curvature effect). The energy evolution follows E ' (t). Using 
F'(t) a more thorough fit of the HIC can be made giving 71, which gives an estimate of 
-R/^sh'/A' or /3 s h'(r s /r) (independent of R). The spectra arise from integrations along 77 = 2 
which, depending on the details of the case, maybe provide a possibility to determine a- in t. 

Case III. Intermediate cases where S-curves are seen (e.g. pulses in Figs. 2a, b, and 
c): The low energy section of the HIC {i.e. at late times) will follow rj = 2 and gives the 
value of r^g and R. With this knowledge the light curve can be deconvolved and F'(t) can 
be found. A fit to the HIC can now be made to find 7Z [which gives Td yn and A'] and 77i nt 
which will be revealed from the early part, [which will allow the determination of E'(t)}. A 
corresponding softening of the spectra as described in the paper should be present. 

BR01 found that in several GRBs, with two separable pulses, the HIC index varied less 
from pulse to pulse in a single burst as compared to its variation in different bursts. This 
requires that the pulses in multi-pulse bursts be produced in shocks created in a similar 
environment, with similar values of R, T TC \, T, 7Z, i] int , n and B. This could happen in a 
scenario in which the two long pulses are created as two similar shells catch up with a leading, 
slower, more bulky shell that has already been significantly decelerated due to interaction 
with the circumburst environment. Such pulses then occur approximately within the same 
environment, at roughly the same distance R (therefore same r ang ) and r re i. This scenario 
also increases the value of r re i, which implies a higher magnetic field, radiative efficiency, 
and a minimum electron Lorentz factor, and a higher synchrotron peak frequency: 

h»s = T^r£ 7c 2 = 5eV(r rcl - I) 2 - 5 1 , (25) 

4nm e c \10 lt> cm/ 

where we have used the relations for B and 7 C described at the beginning of §3. With r re i ~ 
100, his s = 500 keV, so that the expected synchrotron spectrum will peak in the BATSE 
window and will require no additional boost, for instance, from Compton upscattering as in 
the Synchrotron-Self-Compton model (SSC) (Panaitescu & Meszaros 2000). This scenario 
is similar to that of the external shock model normally proposed for the generation of the 
afterglows, which has difficulty to explain the prompt gamma-ray emission because of its 
high variability (Fenimore et al. 1996). However, the GRBs discussed in this paper are 
smooth with few pulses and do not exhibit the high variability of more complex bursts, so 
that this objection is not applicable. 
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